Elastic behavior of spherical nanodroplets in head-on collision 
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Simulation results for head-on collisions of equal-sized spherical polymer nanodroplets using 
molecular dynamics are presented. Elastic behavior of an initial compressed phase for the col- 
liding droplets is analyzed. Deformations and contact radii of the nanodroplets are compared with 
the Hertzian model of elastic solid balls. It is found that at least the initial phase of collision can 
be explained by this continuum model, except at the very moment of the beginning of collision. 
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A nanodroplet has so large a surface-to-volume ra- 
tio that its behavior is quite different from the ordi- 
nary macroscopic droplets which we see in every-day phe- 
nomena. Recently, there are growing interests in nan- 
odroplets collisions since they are closely related 
with applications such as nanofluidics, drug delivery, ink- 
jet printing, etc. In particular, Deming and MasonQ 
made nanodroplets for fighting cancer. Their eventual 
goal is to make possible that the nanodroplets can harm- 
lessly enter into cells and then release medicinal cargo. 
But making better drug-delivery vehicles still needs a 
fundamental understanding of nanodroplets behavior. 



Contact or collision problems in general are very com- 
plicated to model theoretically or even numerically. Here, 
we focus on the head-on collisions of two equal-sized 
spherical polymer nanodroplets using molecular dynam- 
ics. Molecular dynamics (MD) has become an indispens- 
able tool to study the phenomena in nano scale [gj. To 
understand the initial compressed phase of the collision, 
deformations, longitudinal forces and contact radii, etc 
are measured from the simulations. The simulation re- 
sults are analyzed with the Hertzian model j7,]-[8j] of collid- 
ing elastic solid balls. Originally, the Hertz model deals 
with compressed elastic balls such as steel balls. It as- 
sumes Hooke's law between stress and strain. As far as 
I know, there is not yet any attempt to connect nan- 
odroplet collision with the Hertzian theory of elasticity. 
Thus this is a first report to explain the initial phase of 
droplet collision. 



We will first present the Hertzian theory of elasticity 
for equal-sized spherical elastic balls. Next, we describe 
MD simulation method and results. Then, these simu- 
lation results are compared with the theoretical predic- 
tions. Finally, we will summarize our findings. 



When two elastic balls are colliding with each other, 
conservation laws for energy, linear momentum and angu- 
lar momentum are generally applied before and after the 
collision. In this case, the collision process is treated like 
a black box: we ignore the collision process itself. How- 
ever, we want to know what is going on during the colli- 
sion process. This collision process seems to be similar to 
that of the statical compression of two balls, which can 
be explained by the Hertz model. Here, we will summa- 
rize the Hertzian modelQ-Q of contacting elastic balls 
which have an equal-sized spherical shape with a radius 
R. The contact is made by pushing them together. Fig. 
[T^) shows the very moment of contact. Let us denote 
the coordinate system x and x' as being positive in ei- 
ther direction from their contact point. Fig. [Hd) shows 
their squeezed state by a force F„ where F„ is a normal 
component of an applied force from the other ball. The 
deformed surfaces are thus described by equations of the 
form X — x{y,z) for the right ball and x' = x'{y,z) for 
the left ball, respectively, as in Fig. [TJ Let u and u' 
be the x-components of the displacement vectors of the 
contact point on the deformed surfaces of the two balls, 
respectively. They are symmetrically deformed about an 
axis connecting the centers of balls. Then, from Fig. [Hd), 
we have the equality, 



u' = i, 



(1) 



where ^ is a total deformation of the colliding balls. Thus, 
the total deformation ^ varies with the force F„. 

The material is assumed to be isotropic and homoge- 
neous. Linear elastic Hooke's law holds between stress 
and strain. It is further assumed that there is no ad- 
hesion between balls. Furthermore, we consider only a 
small ^ limit. In this limit, the total deformation is 
given by 



(2) 
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where 7 = 3(1 — v'^)/2Y , v is a. Poisson's ratio and F is a 
Young's modulus. The spherical shapes of balls become 
fiat with the collision as shown in Fig. [TJd). The contact 
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a) b) 

FIG. 1; (Color on-line) a) Two equal-sized spherical elastic 
balls with a radius R at the moment of contact, b) A com- 
pressed state. Dotted lines represent the boundaries of the 
two balls at the moment of their contact. They are com- 
pressed with a contact radius a and total deformation 5. The 
balls just begin to touch each other at the points M and M', 
respectively. 

radius a of this flat contact circle is also given by 

a = F^''\^^Y'\ (3) 

III. SIMULATION RESULTS 

Now, we simulate head-on collisions of equal-sized 
spherical polymer nanodroplets. To keep the spherical 
shape of the droplets, we choose polymer nanodroplets [1], 
since a polymer nanodroplet can have less or no evapo- 
ration during the simulation. 

First, let us briefly describe the simulation method. 
The molecular interactions are taken as follows. The 
polymer chains are modeled by a rather abstract but 
well-studied bead-spring model TO]. The chain length of 
a polymer is L = 10. All monomers in the same droplet 
interact with each other through the Lennard- Jones (LJ) 
potential, given by 

VLAr)=Ae[{r/a)-^^-{r/a)-% (4) 

Neighboring monomers in the same chain in addition in- 
teract with the finite extension non-linear elastic (FENE) 
potential, 

VFENE{r) = -^r^log[l - (r/ro)^], (5) 

where kp is a spring constant and ro is a maximum length 
within which the chain can be maintained. In our sim- 
ulations, we choose ro = 1.5 and kp = 30.0. Between 
molecules in different droplets, only repulsive term of LJ 
potential is applied so that the droplets repel each other 
to match with the Hertz model. 



We performed MD simulations by using LAMMPSfT^ 
package with some modifications or additions, if neces- 
sary. Visualizations of the simulation system were done 
by using VMD[121]. Hereafter all the quantities will be 
expressed in the LJ reduced units. 

We prepare two colliding spherical droplets as follows. 
First, we make a polymer melt with chain length of 
L = 10 such that their positions and velocities corre- 
spond to the density of p = 0.5 and temperature of 
T = 0.5. To make a spherical droplet, we connect each 
molecule with a virtual harmonic spring with a small 
spring constant. The other end of each spring is con- 
nected to the center of the simulation box to compress 
the droplet towards its center. Then we can have a spher- 
ical droplet with radius R = 22.87 with a thermodynamic 
state of p = 0.83 and temperature of T = 0.5. After 
the shape of droplet becomes spherical, we turn off the 
connection of harmonic spring. The number of molecules 
composing one droplet is n = 40, 000 so the total number 
of molecules in the system is = 2n = 80, 000. Another 
colliding droplet partner is made with the configuration 
at different time in equilibration process. Thus they are 
at the same thermodynamic state but with different con- 
figurations. Next, they are displaced an appropriate dis- 
tance . Then the two droplets are approaching each other 
with the same speed v to collide, so that their relative ve- 
locity is gn — iv. We choose a significantly larger value of 
cutoff length Vc = 5.0 for LJ interaction, compared to the 
usual value of rc = 2.5, since, otherwise, they will expe- 
rience unwanted forces, when they enter the interaction 
cutoff region from the far distance. 

Since the molecules in different droplets repel each 
other, the droplets are eventually recoiled. We can mea- 
sure their respective recoil velocities and thus the resti- 
tution coefficient e is determined. The coefficient of (nor- 
mal) restitution is given by 



where gn and g'^^ are normal components of relative ve- 
locities of the colliding objects before and after collision, 
respectively. Anyway, e is related with the forces Fn , the 
material properties and the impact velocity etc. The 
simulation results are summarized in Table HI Notice that 
a coefficient of restitution e is strongly dependent on the 
approaching velocity w, since colliding objects are liquid 
droplets, not the solid clusters. The larger the approach- 
ing velocity v is, the smaller a coefficient of restitution e 
becomes in conformance with the earlier results [2| . Their 
recoil velocities reduce significantly, since their energy is 
dissipated in the collision. Thus their coefficient of resti- 
tution e is quite small compared to the solid clusters. 

Now, let us describe the overall collision processes. The 
collision snapshots are shown in Fig. /refoverall at four 
different times at t = 0, 5.0, 27.5, and 75.0. Of course, 
the shapes of colliding droplets significantly change with 
time as expected. At t = 0.0, one droplet has roughly 
a spherical shape with a radius R w 22.09 and the sep- 
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TABLE I: Simulation data of coefficients of restitution e with 
impact velocity v. 



V 


0.1 


0.2 


0.5 


1.0 


1.5 


2.0 


e 


0.437 


0.349 


0.259 


0.174 


0.131 


0.105 



aration gap between two droplets is set at Axo = 5.0. 
They are approaching each other with a relative velocity 
2uo = 3.0. The spheres are unilaterally compressed along 
the longitudinal direction up to i = 6.5, as shown in Fig. 
/refoverallb) . Note that the colliding region is squeezed 
into a flat surface, but the rest still keeps a spherical 
shape. Then the droplet is elongated along the trans- 
verse directions (in y-z plane) up to t = 30.0 as shown 
in Fig. /refoverallc). At around t — 32.0, its shape is 
maximally deformed; it becomes roughly a prolate ellip- 
soid with a length of minor axis ~ 18.15 and lengths 
of major axis Ly Ki 69.57 and ~ 71.22. Then follows 
a recovery returning to its starting shape up to t = 55.0. 
At t = 55.0, the droplets are finally separated and be- 
gin to recoil each other, i.e. they are eventually receding 
with a recoil velocity v'^ w 0.14 with some small acquired 
transverse velocities. Therefore, they are moving at an 
odd angle from each other, even though the collision is 
assumed to be head-on. During the recoil, they may tend 
to return to a spherical shape due to their surface ten- 
sion. In other words, they are in the stress relaxation 
with a long relaxation time. 

Next, to examine the elastic behavior of the collision 
process, we focus on the first part of initial compressed 
phase. To see more detailed motion, we halved the size 
of timestep during initial compressed phase, since the 
particles involved during the collision are moving rather 
rapidly in this phase. In Fig. 21 we present the force 
Fn vs. deformation ^ in the initial compressed phase up 
to t = 5.0. The deformation ^ can be directly measured 
with the standard image analysis. The dotted points are 
the simulation results. The solid line is fitted to the Eq. 

This may come from two reasons: first, the very 
initial contact tip is not a continuous media. They are 
composed of only a small number of molecules. Second, 
the shape of the tip is not spherical, but irregular and 
rough. The fitting is rather good except at small values 
of Fn. From the fitting to the data using Eq. ([2|), we 
obtain a value of 



(27Vi?)i/3 ^ 0.0253. 



(7) 



Finally, to further confirm the validity of the theoret- 
ical prediction for the elastic theory of the initial com- 
pressed phase of a droplet, the contact radius a in Eq. 
([3]) is calculated. The consequent results of the contact 
radius a vs. force Fn are depicted in Fig. [S] The values 
are also measured with the standard image analysis. In 
the same time range, the deformation ^ vs. force Fn in 
Eq. m show the similar behavior as the prediction of Eq. 
([3]). In Fig. [SJ when the radius a is smaller than 6.0, the 
force Fn is shapely increasing. Or, at around Fn = 800, 
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FIG. 2: (Color on-line) Overall longitudinal force vs. time 
t. The droplets are approaching to each other with a velocity 
uo = 1.5. Note the second peak around at t = 11.0. 





b) 




d) 



FIG. 3: (Color on-line) Snapshots of collision progresses at 
different times, a) t = 0, separated with gap of Axq ~ 5.0 be- 
tween two leading edges, approaching with a impact velocity 
2uo — 3.0, b) t = 5.0, initially compressing stage, still main- 
taining a spherical shape, c) t = 27.5, transversely expanding, 
prolate ellipsoidal shape, d) t = 75.0, receding stage, in the 
process of returning to its spherical shape. 



a sharp increase of stretched S-shape is observed. But at 
around a — 10.0, they show a linear increase. From the 
fitting to the data using Eq. ([S]), we can again obtain a 
value of 



(7i?/2) 



1/3 



0.679. 



(8) 



From the Eqs. ([7]) and ([S]), we can get the estimated 
value of Rest ~ 36.45, which is rather larger than our 
initially setting value of i? = 22.87, but the order of 
magnitude is comparable. 
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FIG. 4: (Color on-line) Deformations ^ vs. longitudinal force 
Fn in the initial compressed phase of a collision. Solid line is 
a fit to the data by Eq. ^ with = 0.73 and (7^2//?)^/^ = 
0.0253. Note the poor agreement near the origin. 
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FIG. 5: (Color on-line) Contact radius a vs. longitudinal 
force Fn in the initial compressed phase of a collision. Solid 
line is a fit to the data by Eq. ((3| with — 0.78 and 
(7-R/2)^/'^ = 0.679. Note again the poor agreement near the 
origin. 
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